Submitted to the IEEE TRANSACTIONS ON INFORMATION Theory, April 2010 



The Re-Encoding Transformation in Algebraic 
List-Decoding of Reed-Solomon Codes 

Ralf Koetter, Fellow, IEEE, Jun Ma, Senior Member, IEEE, and Alexander Vardy, Fellow, IEEE 



O 



m 






> 

m 

in 

vn 

o 
o 



X 



Abstract — The main computational steps in algebraic soft- 
decoding, as well as Sudan-type list-decoding, of Reed-Solomon 
codes are bivariate polynomial interpolation and factorization. 
We introduce a computational technique, based upon re-encoding 
and coordinate transformation, that significantly reduces the com- 
plexity of the bivariate interpolation procedure. This re-encoding 
and coordinate transformation converts the original interpolation 
problem into another reduced interpolation problem, which is or- 
ders of magnitude smaller than the original one. A rigorous proof 
is presented to show that the two interpolation problems are in- 
deed equivalent. An efficient factorization procedure that applies 
directly to the reduced interpolation problem is also given. 



I. Introduction 

REED-Solomon codes are still widely used in digital com- 
munications and data storage. Conventional hard-decision 
Reed-Solomon decoders correct up to {n — k) /2 symbol errors 
for an {n, k) Reed-Solomon code. During the past decade, sev- 
eral breakthroughs in improving the error-correction capabil- 
ity of Reed-Solomon decoders have been achieved. First, Su- 
dan [20] proved that list-decoding of Reed-Solomon codes can 
be viewed as a bivariate interpolation problem, thereby correct- 
ing more errors than previously thought possible. Specifically, 
for an {n,k) Reed-Solomon code, Sudan's algorithm [20] pro- 
duces all codewords whose distance to the received hard-deci- 
sion vector is at most n — \/2kn. The second major step was 
taken by Guruswami and Sudan [8]. Their work showed that 
one can correct even more errors by interpolating through each 
point not once, but m times, where m is an integer parameter. 
For m ^- 00, the list-decoding algorithm of [8] corrects up to 
n — ^/nk errors. These results were later extended by Koet- 
ter and Vardy [12], who showed how the interpolation multi- 
plicities in the algorithm of [8] should be chosen to achieve 
algebraic soft-decision decoding of Reed-Solomon codes. 

Both list-decoding and algebraic soft-decision decoding rely 
upon interpolation and factorization of bivariate polynomials. 
Such interpolation and factorization is much more computa- 
tionally intensive than the Berlekamp-Massey algorithm used 
in conventional hard-decision decoding. Various efficient algo- 
rithms for bivariate interpolation and factorization have been 
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proposed by Alekhnovich [2], Augot- Pecquet [3], Feng-Gira- 
ud [4], Gao-Shokrollahi [5], Koetter [9], Lee-O' Sullivan [13], 
Nielsen-H0holdt [16], Olshevsky-Shoki-ollahi [17], Roth-Ruck- 
enstein [18], Wu-Siegel [24], among others. While polynomial- 
time, these algorithms fall short of making the required com- 
putation feasible in practical applications that usually involve 
long high-rate Reed-Solomon codes. 

In this paper, we present a series of transformations that 
drastically reduce the space and time complexity of the bivari- 
ate interpolation process, by a factor of at least n^ /{n—k)-^. 
These transformations were first reported in the conference pa- 
pers [6], [10], and [11]. Later, the resulting re-encoding and co- 
ordinate transformation procedures were used, without proof, 
in [1,7,14] and other papers. Our main goal herein is to pro- 
vide a streamlined formulation of this transformation process, 
and of the corresponding factorization procedure, including ri- 
gorous proofs and detailed examples. 

The rest of this paper is organized as follows. In Sectionllll 
we introduce the necessary definitions and establish (in Theo- 
rem[T]l a basic property of the birational isomorphism that un- 
derlies the coordinate transformation process. The bivariate in- 
terpolation problem encountered in list-decoding and in alge- 
braic soft-decoding of Reed-Solomon codes is briefly reviewed 
in SectionUn] We also review in Section|lII] the Grobner-basis 
bivariate interpolation algorithm due to Koetter [9]. The re- 
encoding and coordinate transformation complexity-reduction 
method is described in detail in Section|lV] This section also 
contains a rigorous proof of the validity of this method. In Sec- 
tion|V] we show how the last decoding step, namely bivariate 
polynomial factorization, should be carried out in order to take 
full advantage of the reduced interpolation procedure of Sec- 
tion|lV] We conclude with a brief summary in SectionlVTl 



II. Definitions and Preliminaries 

Let F^ be the finite field with q elements. The ring of polyno- 
mials over Fq is denoted F^j [X] . Reed-Solomon codes are ob- 
tained by evaluating certain polynomials in F^ [X] on a set of 
points V C F^. Specifically, the Reed-Solomon code Cq{n,k) 
of length n and dimension k is defined as follows: 

Cq{n,k) t^ {(/(^i) /(^„)) : ^1,^2 ^neV, 

/(X)eF,[X], deg/(X)<fc} 

The set D is usually taken as F^ or as F * , where F * is the set 
of all the nonzero elements of Wq. Unless stated otherwise, we 
shall assume that T> = F* so that n = q — 1. As in [8,12,18], 
we define the weighted degree of a polynomial as follows. 
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Definition l.Let^(X,Y) = YT=oLjLoai,jX'YJ beabivariate 
polynomial and let ivx, Wy be arbitrary real numbers. Then the 
{zvx, iVY)-weighted degree of A{X, Y) is defined as the maxi- 
mum over all real numbers izvx + jwy such that a, ,■ ^ 0. 

For reasons that will become clear later, we do not restrict 
the definition of weighted degree to the usual case [8,12,18] 
where ivx, Wy are nonnegative integers. Thus the weighted 
degree of a polynomial A{X,Y) can assume negative values. 

Let K^ a denote the ring of rational functions in Wq{X,Y) 
without poles at the point (a, /5) G F^ x F^. A rational function 
A{X,Y) E K/^R always has a power-series expansion in basis 
functions of type (X — a)'(Y — /5)^. Thus we can write 



Ax,Y) = ^^fl,.^.(x-^)'(y-/3y (i) 

i=0;=0 



and define 



coef(^(X + a,y + ^);X'Y^) 



def 



"■',! 



(2) 



Definition 2. The function A{X, Y) is said topass through the 

point {oL, [i) with multiplicity m ifciji = for all i + j < m in 
the power-series expansion ([T])- Further, the multiplicity func- 
tion fia,B '■ ^a,e, ~^ IN is defined as follows: 

lia,i}{A{X,Y)) = max{;« e N : Oj^j = Vf + ; < m} 

The following properties of the multiplicity function follow 
from the well-known fact that /;„, « is an (exponential) valua- 
tion. Given any A{X,Y) and 13{X,Y) in -K^, r, we have 

li^,li{A{X,Y)) = ^o,oi■^{X + a,Y + ^)) (3) 

H^,li{AB) = ii^,fi{A) + }i^,p{B) (4) 

}i,,piA + B) ^ mm{ii,^i,iA),^,,f,{B)} (5) 
Assuming that ^(X, Y)/6(X,y) E K^^p, we also have 

F^,li [sj ^ HocfiiA) - }i^,ii{B) (6) 

Definition 3. Fix a polynomial g{X) over Fg, and let 3 be the 

set of all a E F^ such that g{a) = 0. We define the following 
pair of mappings between points in (F^— 3) x F^; 



f- i^'V) ^ (x, 



y 



g{^) 



cp^: (x,y) H^ {x,ygix)) 



(7) 

(8) 



It is easy to see that (pg and cp„ ^ are birational isomorphisms, 
and (p~^ is indeed the inverse of q>n. 

Definition 4. Fix a polynomial g{X) over Fq. Then the maps 
0g and 'Pg^ from F^ (X, Y) to itself are defined as follows: 



OO OO CO Y^OO -y-/ 

j=oi=o j=o S[^)' 



(9) 



Thus <I>g, respectively <I>g \ maps a rational function A{X,Y) 
into A{X,Yg{X)), respectively A{X,Y/g{X)). 



The following theorem establishes an important property of 
the map <J'^ that will be crucial for our results in Section llVl 

Tlieoreml. Let (a, /S) E (F,-3) x F,,andJet^(X, Y) eX^^p. 
Define (a, 7) = (pg{oi, ^) and B{X,Y) = <Pg{A{X,Y)).Then 

li^,fi{A{X,Y)) = f/„,^(B(X,Y)) 

Proof. Since, by assumption, A{X, Y) does not have a pole 
at (a, /5), and6(a, 7) = ^(a, 7g(a)) = ^(a, j6), we see that 
B{X, Y) does not have a pole at (a, 7). Therefore, we can ex- 
pand B{X, Y) in basis functions (X — a.y{Y — 7)-' as follows: 

B{x, y) = E E hi (X - ^y{y - ly m 

i=Oj=0 

With the coefficients bjj determined by ( fTOt . our definition of 
the multiplicity function implies that 

li,,^{B{X,Y)) = min If + ;■} 

Since the maps <I>g and (Pg^ in (|9|l are inverses of each other, 
we can express A{X,Y) as follows: 



A{X,Y) = 0g'{B{X,Y)) = b(x, 



Y 



g{X) 



(11) 
(12) 



i=0;=0 



gix) J 



Dividing X — a into ^(X) with remainder g{ix.), we can write 
^(X) = (X-a)^(X)+^(a)and 

7g(X + a) = 7X/z(X + a)+7g(fl:) = li + jXh{X + a) 

for a nonzero polynomial h{X). With this, it follows from (O, 
Q, and ([13]) that 



Fa,/3(-4) ^ min <^ fi^^p (X - a)' 



.YY- 



mm I ^0,0 (x' 
mmLoJx' 



lg{X) \'" 



g{X) J 
Y+^-jg{X+a)\' 



giX+^) 
Y-jXh{X + oc)'^' 

g{X + c^) 



Since g{a) ^ 0, we can use dU and (|6]l to further simplify the 
right-hand side above as follows: 



f 0,0 X 



Y-jXh{X + oc)\' 



g{X + a) J J 

z>o,o(X) + m,o (y - jXhiX + a)) - j^o,o {g{X + a) 

It is clear that }iofi{X) = ?^io,o(y - jXh{X + a)) = 1, while 
Fo,o{g{X + d)) = 0. Hence, we finally conclude that 

^,,|,{A{X,Y)) ^ miii{z+y} = fi,,^(fi(X,Y)) (14) 

To establish the inequality in the other direction a similar ar- 
gument can be used. First expand A(X, Y) in basis functions 
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of type {X — a)' {Y — fi)J in order to define the coefficients Oj; 
such that //„ fi(^(X, y)) = mirifl, ^^o{i +/}■ Then write 

B{X,Y) = 0{A{X,Y)) = A{X,Yg{X)) 



J2J:a,j{X-ay{Yg{X) 

i=Oj=0 



to conclude that 



{?/o,o(x'((Y + 7k(X + ^)-/5y)} 



min 



But(Y + 7)g(X + a)-^ = yg(X + a)+7X/2(X + a)andit 
is easy to see that ^o,o(Yg(X + a) +7X/!(X + a)) = 1. This 
shows that ^fl.^^(S(X, Y)) ^ mirifl. .^o{^+/} ^"d, in conjun- 
ction with (fl4l i. completes the proof of the theorem. | 



III. The Interpolation Problem 

Our interest in the foregoing definitions and results is mo- 
tivated by the fact that, as a consequence of Bezout's the- 
orem, two polynomials ^(X, Y) and B(X, Y) cannot both 
pass with high multiplicity through an arbitrary large number 
of points without having a common factor. In particular, the 
polynomial Y — /(X), with deg/(X) < fc, passes through the 
n points (^i,ci), (^2,C2), • • ., (^«,c„), where C; =/(?,) may 
be thought of as the n transmitted symbols. Then Bezout's the- 
orem implies that any nonzero polynomial Q,(X,Y') such that 

Y^n-AQ) > degi^,_i Q(X,Y) 

is divisible by Y — /(X). This observation leads to the interpo- 
lation-based decoding algorithms of [8,20] and [12]. The cen- 
tral idea of these decoding algorithms is to construct a polyno- 
mial Q{X,Y) that passes through a prescribed set of points 
def 



V 



^{xi,yi),{x2,y2), ■ ■ ■ , {xs,ys)} 



(15) 



where Xi, X2, . . . , x, G X* and 1/1,1/2/ • ■ • / J/s G Fg, with the pre- 
scribed multiplicities wixi,i/i/ wz 1-2,1/2' • ■ • ' ^x,,y,, which are pos- 
itive integers. If these points and multiplicities agree "suffi- 
ciently well" with the n points (^/, C/) that determine the trans- 
mitted codeword, then divisibility of Q{X, Y) by Y — /(X) is 
guaranteed. The reader is referred to [12] for more details. 

In all cases, a key part of the decoding algorithm consists 
of solving the following interpolation problem. 

Definition 5 (Interpolation problem). Given a set of s points 
V = {{xi,yi), (x2,i/2)/- ■ •/ (^s,i/s)} and a set of s multiplici- 
ties M = { WJ-Ti ,1/1 / wz 12,1/2 ' ■ • ■ ' ^XsAjs}' f^s interpolation prob- 
lem consists of computing a polynomial Q{X,Y) ^ such that 

}ix„y,{Q{X,Y)) ^ m^.,y^ for all (x„y,) e V (16) 

and degj i^_iQ{X,Y) is minimal among all bivariate polynom- 
ials that satisfy the interpolation constraints ( fT6] l. 

We shall refer to this interpolation problem as IP^ j.-i {'P/ M), 
and say that Q{X,Y) is a solution to IPi j^_i{V, M) . Observe 
that Xi and X; in the set V = {{xi,yi), {x2,y2), ■ ■ ■ , (^s,i/s)} 
do not have to be distinct, all we require is that they belong 



to T>. In fact, in soft-decision decoding, we often interpolate 
through different points having the same X-coordinate [12]. 

By definition (see Definition|2]i, requiring that a polynomial 
Q(X, Y) passes through a given point with multiplicity m im- 
poses ^l^m{m+V) linear constraints on the vector space of 
polynomials in two variables. Hence, solving the interpolation 
problem IPj fc-i ("P, M) is tantamount to solving a system of 



N(M) = ^Emx.y,^,,^,^!) 



(17) 



1=1 



homogeneous linear equations. As shown in [8], [12], and other 
papers, there are precisely 



,^1,/c- 



-IW = 






k-1 



k-1 



(18) 



monomials X'Y^ with i + {k—l)j ^ S. Hence, choosing S to 
be large enough will guarantee a solution to IPj j.-i(^/M). 
Let S* be the least integer such that Xl,k-l{^*) > N(M). Then 
the weighted-degree of a solution Q{X,Y) to IPij-_i(P, M) 
satisfies deg^ f-_^ Q{X, Y) ^ S*, and its Y-degree is at most 

3* 



def 



k-1 



(19) 



In principle, IP^ if_i{V,M) is a simple linear problem that 
can be solved in a number of ways. Numerous algorithms for 
this purpose have been proposed in [1,2,4,6,7,9,13,16-18,22]. 
In what follows, we briefly review the interpolation algorithm 
of Koetter [9], which is widely recognized as one of the most 
suitable for implementation in practice. VLSI architecture for 
this algorithm has been developed in [1,6,7] and other papers. 

In order to describe Koetter's interpolation algorithm [9], we 
first need to extend the notion of (1, A:— 1) -weighted degree to 
a monomial order. Explicitly, we say that X''Y ~<k X'Yi iff 

a + {k-l)b < i + {k-l)i 

or (20) 

a + {k-l)b = i + (fc-1); and b < j 

Every polynomial in F^j [X, Y] now has a well-defined leading 
term under ^fc, and we can impose a total pre-order ^/t on po- 
lynomials in F(j[X, Y] by comparing their leading terms. 

Notably, Koetter's algorithm computes much more than just 
a solution to the interpolation problem IPj j--i(P, M). It pro- 
duces a set of polynomials, which forms a Grobner basis for 
the ideal of all polynomials over Fg that satisfy the interpola- 
tion constraints defined by V and M. This Grobner basis ^ 
consists of r + 1 polynomials: ^ = {Qq, Qi, . . . , Qr}, where r 
is given by ( fT9] l. The computation is initialized by setting 

^"" := {1,Y,y2,...,Y'} (21) 

Upon initialization, Koetter's algorithm goes through N{M) it- 
erations, imposing each of the N{M) linear constraints in ([TtT i 
one-by-one. For example, suppose that at iteration ; of the al- 
gorithm, we are deaUng with the constraint 

coef(Q(X + x,Y + y);X''Y*') = 

for some {x,y) E V and some a + b < nix^y- Then, given the set 
^ ~ { ^0 ' ^1 ' ■ • ■ ' '3r } computed after the first i — 1 
iterations, the f-th iteration consists of the following. 
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Update Basis (^; {x,y);a,l 



1 1 1 For all i = 0,1, . . . ,r, compute the discrepancy of Q:' 
which is given by: 



def 



ay/h\ 



Aj = coei{g-' '{X + x,Y + y);X''r 

If Aj = for all j = Q,l,...,r, set ^'" := ^'"" and stop. 

[2] Among t/Q^",^!'"", . . .,t/ji'"", find the least with resp- 
ect to -<fc polynomial such that its discrepancy is nonzero. 
Let Qf' (X, y) denote this pivot polynomial, so At 7^ 0. 

|3 I For all j = 0,1, . . .,r, except ;' = t, compute 

g['ix,Y) := gl'\x,Y) - ^ gr''(x,Y) 

Then update the pivot polynomial Qf ' (X, Y), namely set 

g';\x,Y) := (x-x)g;'-''(x,Y) 

It can be shown [22] that for all i, the Y-degree of the lead- 
ing monomial of Q- (X, Y) is exactly /. Along with ( l20b . this 
implies that the polynomials Q^ , Q^ , . . . ,Qy always have dis- 
tinct orders with respect to ^jt and guarantees that the choice 
of the pivot polynomial in Step [2] is unique. 

Using the Update Basis (^; {x, y);a, b) procedure above, 
the entire interpolation algorithm can be formulated as follows. 

Koetter's Interpolation Algorithm 

Input: A set of points V = {{xi,yi),{x2,y2), ■ ■ ■ ,{^s,ys)} 
along with their multiplicities mxj^,yj^, nix2,y2' ■ ■ ■ ' ^Xs,ys- 

Initialization: Set ^ := { 1, Y, Y^, . . . , Y"" } . 

Iterations: For all {x,y) e V, do the following: for a := to 
mx,y — 1, then for b:=0 to nix^y — a — 1, update ^ using 
the procedure Update Basis ($^; {x,y);a,by 

Output: Return the least with respect to -<t polynomial in ^. 

It is not difficult to see that the number of additions and mul- 
tiplications in F^ required to solve IPi;c_i(P, M) using this 
algorithm is 0{rN'^), where N = N{M) is the total number 
of linear constraints given by ( fTTI i. While this is substantially 
faster than straightforward Gaussian elimination, the problem 
is that the number of equations N is often too large to make 
an 0{rN-^) computation feasible in practice. The following ex- 
ample sheds some light on the magnitude of this problem. 

Example la. Let Cq{n,k) be a Reed-Solomon code of length 
n = 255 and dimension k = 239 over F256. A typical interpo- 
lation problem arising in algebraic soft-decision decoding [12] 
of Ccf{n,k) might involve the following multiplicities: 



(22) 



for a total of 6192 linear equations. Using ( fTSl ). we find that the 
corresponding value of S* is 1598. The Y-degree of a solution 
Q{X,Y) to this interpolation problem is r = 6, in view of ( fT9] l. 



multiplicity 


7 6 5 4 3 


2 1 


# of points 


229 12 10 4 3 


10 10 


# of constraints 


6412 252 150 40 18 


30 10 



Estimating the complexity of Koetter's interpolation algorithm 
as rN^, we arrive at roughly 230 xlO^ operations. In fact, the 
algorithm requires exactly 159.56 xlO^ finite-field multiplica- 
tions. This figure is precise: it was obtained by actually imple- 
menting the algorithm, and counting the number of finite-field 
multiplications in software. The fastest algorithm we found for 
the particular interpolation problem in (|22]) is due to Lee and 
O' Sullivan [13]. This algorithm takes exactly 45.37x10^ finite 
field multiplications, which is still prohibitive. □ 

This example illustrates a major problem with interpolation- 
based decoding. While, for a fixed maximal multiplicity, the 
decoding complexity is bounded by a polynomial in the length 
of the code, the actual complexity of solving the interpolation 
problem is prohibitively large in practice. In the next section, 
we present an algorithmic transformation that drastically redu- 
ces this complexity. Before describing this transformation in 
detail, let us use the following example to illustrate the savings 
in complexity that can be achieved with this method. 

Example lb. Consider again the interpolation problem in (l22l i. 
Judiciously choosing the re-encoding point set, we can elimi- 
nate the k = 239 points with the highest multiplicities: the 229 
points of multiplicity 7 as well as 10 of the 12 points of multi- 
plicity 6. This leaves only 290 linear equations to solve, rather 
than the original 6912. This is a much more feasible task. In 
fact, using Koetter's algorithm to solve the reduced interpola- 
tion problem requires only 350x10^ finite-field multiplicati- 
ons. The reduction in complexity, by a factor of 456, is aug- 
mented by a corresponding reduction in memory requirements, 
due to the fact that the polynomials operated upon during the 
interpolation procedure have much small degree. Specifically, 
the largest X-degree we need is about 40, instead of 1598. □ 



IV. A Complexity Reducing Transformation 

Rather than seeking an efficient way to solve IPi j--i(^/M), 
we will modify the interpolation problem itself, by means of 
a shift and a coordinate transformation. Our approach is sim- 
ilar to the re-encoding idea of Berlekamp and Welch [23]. 

A. Re-Encoding and Shift 

Given the point set V = {(xi,i/i), (x2, 1/2)/ • • •/ (^s/i/s)}, we 
first identify some k points (x,^,:/,^), (x,2,i/;J, . . ., (x,-^,i/,J 
in V such that Xj^,Xi^, . . .,Xi E V are all distinct. Define 

def 



7^ 



{{Xi,,yn),{xi2,yi2),...,{xi^,yi,)} 



(23) 



Observe that if the set V contains n — e < k points with dist- 
inct X-coordinates, then the resulting interpolation polynomial 
Q(X, Y) will have at least q'-^"-^) factors of type Y-/(X). 
This situation corresponds to e > n — k erasures, in which case 
the transmitted codeword cannot be uniquely determined. This 
shows that, unless the interpolation problem IPj j._i(P, M) is 
ill-conditioned by too many erasures, a re-encoding point set 
TZ with the required property always exists. 

Note that there will usually be exponentially many ways to 
choose TZ from V. As far as the theory developed in this paper 
is concerned, the choice of TZ is arbitrary. In practice, the set TZ 
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will be chosen to consist of the points with the highest possi- 
ble multiplicities (cf. Example 1). To simplify notation in what 
follows, we assume without loss of generality that TZ consists 
of the first k points of V, that is 



n 



|(xi,1/l),(x2,1/2),...,(Xfc,l/fc)| 



The re-encoding point set TZ determines the unique re-encoding 
polynomial e(X) of degree < k, defined by 



e{xi) = yi for all {xi,yi) E TZ 



(24) 



Observe that the codeword c' obtained by evaluating e(X) at 
^1, f 2/ ■ • • / ?n agrees with the "given" values i/i, 1/2, ... , yt at 
the k positions corresponding to Xi, X2, ■ ■ ■ , X/^. Thus comput- 
ing e(X) is equivalent to re-encoding through k given values 
at some k positions. If these k positions are consecutive and 
Cq{n,k) is cyclic, this can be achieved through division by 
the generator polynomial for Cq{n,k). Otherwise, such re-en- 
coding is tantamount to correcting n—k erasures in Cq{n,k). 
Various efficient algorithms for this purpose are known [15,23]. 
Given the set V = {{xi,yi), {x2,y2),- ■ ■ , (^s/i/s)} and the 
re-encoding polynomial e(X), we define 



V 



' = |(xi,i/i-e(xi)),(x2,i/2-e(x2)),...,(xs,i/s-e(xs))| 

Notice that, by the definition of e(X) in ( |24] |. the first k points 
in V' are of the form (xi,0), (x2,0), . . ., (x/f,0). 

Lemma 2. Consider an arbitrary point (a, fi) eV and the cor- 
responding point (a, /5') eV, where fi' = /5 — e(a). Furttier, 
let A{X, Y) and B{X,Y) be arbitrary polynomials related by 

B{X,Y) = ^(X,y + e(X)) (25) 

A{X,Y) = B{X,Y-e{X)) (26) 

li,,l,{A{X,Y)) = y,^p{B{X,Y)) (27) 

Proof. The proof is very similar to the proof of Theorem[T] 
Expand B{X,Y) in the basis (X - a)'(y - ^'y to write 

00 00 

B{X,Y) = E E hi (X - ^nr - fi'y (28) 

1=0 ;=0 



Then 



so that j,inff{B{X,Y)) = min;,, .-^oi'+i}- Then use 
and the fact that e(X) = (X — oc)h{X) + e{a) for some poly- 
nomial h{X), to express A{X + a,Y + fi) as follows: 

A{X + cc,Y + ^) = B{X + cc,Y + ^-e{X + oc)) 
i=Oj=0 

00 00 ; 

= ELhi^'[y-XHX + ci)) (29) 
i=Oj=0 

Since ^0,o(X) = piofi{Y — Xh{X+oc)) = 1, it follows immed- 
iately from ©, (gji, Q, and (|25]l that 

li,,l,{AiX,Y)) ^ miii{/ + /} = y,^p\B{X,Y)) 

The proof of the inequahty y^^f^, {B{X, Y)) ^ y^^fi {A{X, Y)) 
is essentially identical to the above, and is omitted. | 



Remark. It is not a coincidence that the proofs of Theo- 
rem[T] and Lemma|2] are so similar. Both results can be estab- 
lished in the general framework of birational maps between al- 
gebraic varieties [19]. However, there is no real need to invoke 
the heavy machinery of birational isomorphisms in this paper 

Theorem 3. A polynomial Q{X,Y) is a solution to the inter- 
polation problem IPx yt-i (^/ A/T) if and only if the polynomial 

Q'{X,Y) = Q(X,Y'+e(X)) is a solution to lP^^^__iiV',M). 

Proof. Let J(P, M) denote the ideal of F^ [X, Y] consisting 
of all polynomials A{X, Y) that satisfy the interpolation con- 
straints determined by V and M, namely such that 

y:c,y{AiX,Y)) ^ m:c,y for all {x,y)EV (30) 

It follows from Lemma|2]that Q(X, Y) E T{V, M) if and only 
if S'(X,Y) = Q(X,Y + e(X)) is in J(P',M). It remains to 
prove that Q(X, Y) is of minimal (1, fc—l) -weighted degree in 
its ideal if and only if so is Q'{X,Y). Observe that 

degi^,_i Q'(X,Y) = degi^,_i Q(X,Y) 

This follows from the fact that dege(X) ^ fc — 1, and there- 
fore deg^ )c-l^ ~ ^^Zi /c-i (^ + ^(■^)) = k — 1. In general, it 
is easy to see that the transformations in ( |25] l and ( |26] l preserve 
the (1,A:— l)-weighted degree. Now, assume to the contrary 
that Q'(X, Y) is minimal in I{V', M), but I{V, M) contains 
a polynomial A{X,Y) such that 

deg,,_,A{X,Y) < degi,_iS(X,Y) 

But then we can use (|25] | to produce a polynomial B{X, Y) in 
X(V' , M) whose (1, A:— l)-weighted degree is strictly less than 
that of Q'(X, Y), a contradiction. | 

It follows from Theorem[3]that instead of solving the interpo- 
lation problem IPj j^_i{V,M) directly, we may first compute 
a solution S'(X,Y) to IPi,fc_i(P',M) and then set 

S(X,Y):=Q'(X,Y-e(X)) 

Every solution Q{X,Y) to IP^ ^_i{V ,M) can be obtained in 
this manner. The next example illustrates the process in detail. 

Example 2a. Let F^ = Fg = {0, 1, a, a^, oc^, oc^, cc' , a^}, where 
a is a root of the primitive polynomial X + X + 1. Take fc = 2 
and T) = {1, a, cp-, 0?}, and consider the Reed-Solomon code: 

C8(4,2) = {(/(l),/W,/(a2),/(a3)) :/(X)=fl + ^ 

Suppose that the codeword c = (1, a^, cc" , a), corresponding to 
/(X) = a^ + a^X, was transmitted. Further assume that, upon 
observing the channel output, the multiplicity assignment algo- 
rithm of [12] produces the following interpolation problem: 

{x,y) (a,a4) (a^,^^) (a^,^^) (a^,!) (a^,^) (l,a) (1,1) 



m 



x,y 



1 



1 



This roughly corresponds to receiving the hard-decision vector 
v_ = (a, a^, c& , 1), with errors in positions 1, 2, 4. Since k = 2, 
the (1, A:— 1) -weighted degree reduces to the total degree in this 
case. Hence, we have ^(^i 1 {S) — {6 + l){5-\- 2) /2. Since there 
are N{M) = 9 linear constraints, the corresponding value of S* 
is 3, and r = 3 by (fT9l l. Koetter's interpolation algorithm thus 
initializes the Grobner basis as ^'"^ = {1, Y, Y^, Y^} and pro- 
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ceeds through the 9 iterations detailed in Table|T] Note that, at 
each iteration, the Grobner basis polynomials are arranged in 
TableUin ascending order with respect to -<i:. Thus the output 
of the algorithm is Q{X,Y) = g^'''(X,y), given by 

(1 + oc^X + uX^) + (a^ + X + X^)Y + (a^ + X)Y^ (31) 

It is easy to verify that Q{X,Y) satisfies all the interpolation 
constraints. Also, the polynomial Y — (a^ + oc^X), which cor- 
responds to the transmitted codeword, is a factor of Q{X,Y). 
In fact, the complete factorization of Q{X,Y) is as follows: 



Q(x,y) = (a3 + x)(y- 



-x))(y- 



^X)) 



Now, let us apply the shift transformation of Theorem[3]with 
respect to the re-encoding point set 

7^ = {{a,ci^),{a.la^)} (32) 

which consists of the first k = 2 points of V. This set TZ de- 
termines the re-encoding polynomial e(X) = a^ + a^X, with 
e(a) = a* and e(a'^) = a^. The polynomial e(X), in turn, de- 
termines the shifted interpolation point set V' given by 

(x,y) (a,0) (a2,0) (a^,^*) (a^,^) (a^l) (1,0) (l,^^) 



m^cy 2 1 1 1111 

The Koetter algorithm for the shifted problem IPx k-i {'P' , M) 
proceeds as shown in in Tablellll and produces the output 

Q'{X, Y) = (a^ + X + X^)Y + (a^ + X)Y^ (33) 

Again, it is easy to verify that Q'{X, Y) satisfies all the interpo- 
lation constraints of IP^ j._i(P', M), and factors as follows; 

S'(x,y) = (a3 + x)y(y-(a + x)) 

We see that indeed, S(X,y) = S'(X,y- (a^ + a''X)) as ex- 
pected. Note that we could also obtain the y-roots of Q(X, y), 
namely Y — {a^+a^X) and Y — (a^ + a^X), by first factoring 
Q'(X, y) and then shifting its y-roots. □ 



B. Coordinate Transformation 

We are now ready to proceed with the complexity reducing tran- 
sformations. For m G Z, define [m]+ = max{m, 0}. 

Lemma 4. The polynomial A{X,Y) = EjLoCj{X)Yi passes 
through a point (a, 0) with multiplicity m if and only if the uni- 
variate polynomials cii{X) are all divisible by (X — a) ''"^■'1 . 

Proof. Expand aj{X) in the basis functions (X — a.)!, that is 
write fl;(X) as aj{X) = J^'fLo^hji-^ ~ '^)'- Then the expans- 
ion ([T]) of ^(X, y) at the point (a, 0) is given by 

^(X,y) = Y^a,,j{X-ocy{Y-Oy = J^ a,,j{X - ayY' 

i,j=0 i,j=0 

Clearly, the polynomial fl;(X) is divisible by (X — a)'™"^! if 
and only if a, , = for all nonnegative i < m — j. This is just 
a reformulation of the definition of multiplicity. | 

For ease of notation, let us henceforth denote the multiplic- 
ities 'W;cj yj, ffij:2,y2/ ■ • • / '^xj-m °f '^he points in the re-encoding 
set 7?. as Vj, 1/2, ... , Vj^. Thus i/,- = m^j^yi for all f = 1, 2, . . . , fc. 
The following corollary is immediate from Lemma|4] 



Corollary 5. The polynomial A{X,Y) = Y^f^^ a j{X)Yi passes 
through the k points {xi, 0), (x2, 0), . . . , (xjt/ 0) with multiplic- 
ities Vi,V2, ■ ■ ■ , V]f if and only if all the polynomials aj(X) are 
divisible byYti=i{X - X;)!^'"^!^. 

From Corollary|5] we conclude that a solution Q'(X, y) to 
the interpolation problem IPj if^i{V' , M) must have the form 

Q'{X,Y) = X: (qi{X)Yl{X-x,)^^>-iAYi (34) 

for some polynomials qQ[X) , qiiX) , . . . ,qr{X). This observa- 
tion makes it possible to reduce the number of iterations in the 
interpolation algorithm by pre-solving for the first k interpo- 
lation points (xi,0), (x2/0), . . ., (x|t, 0). Specifically, in lieu 
of ( l2n i. we initialize the Grobner-basis polynomials as follows: 

gf\X, Y) = y^' n(X - X,) '^' "^'^^ for ;• = 0, 1, . . . , r (35) 
(=1 
By Corollary|5] these polynomials already satisfy all the inter- 
polation constraints determined by (xi,0), {x2,Q), ■ ■ ■ , (xj., 0). 
Thus it remains to enforce only those interpolation constraints 
that correspond to the other points in V' . 

Example 2b. Consider again the re-encoding transformation in 
Examplel2al Recall that k = 2, and the relevant interpolation 
points are (xi,0) = (a, 0) and (x2, 0) = (a^, 0), with multipH- 



(36) 



cities 1/1=2 and 1/2 = 1 respectively. According to (l35l l. let us 
initialize the Grobner basis as follows: 

g'°\X,Y):={X-ocf{X-oi^) 
g';\x,Y) :=y(X-a) 
g!^\x,Y) :=y2 
g'^\x,Y) :=y3 

Comparing with Tablellll we see that these are precisely the po- 
lynomials produced by the Koetter algorithm after the first 4 it- 
erations, starting with the standard initialization {1, y, Y^, Y^}. 
Thus, given the initialization in (|36] |. it remains to proceed with 
the last 5 iterations of the algorithm, exactly as before. n 

This is nice but not enough, since most of the computation 
takes place in the later iterations of the interpolation algorithm 
(cf. Tables I and II). To further reduce the interpolation comp- 
lexity, let us introduce the auxiliary polynomials ^(X), i/'(X), 
and the "tail" polynomials tAX), defined as follows: 

^(X) tf n(X-xO (37) 



^{x)^MyI{x- 

(=1 

tj{X)^^f\{X- 
/=1 



J;-i',]^ 



(38) 



for ;■ = 0, 1, . . . , r (39) 



Combining these definitions with (l34l i. we can express a solu- 
tion Q' (X, y) to the interpolation problem IPj j._i (P', M) as 

S'(X,y) = ^p{X)j^qj{X)tj{X) (^)' (40) 

Computing Q'(X, y) and thereby solving both IPi,jc-l (^'/ A4) 
and IPj ^_i{V, M) (in view of Theorem[3]l thus reduces to find- 
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ing qQ(X),qi{X),. . .,qr{X) in ( |40] i. The following two pro- 
positions show that computing qQ{X), qi (X), . . . , qr{X) is tan- 
tamount to solving a much smaller interpolation problem! 

Proposition 6. Let {a, /5) e V be such thatg{a) 7^ 0. Then the 
polynomial Q' (X, Y), as defined in ( |40] |. passes through (a, j6) 
with multiplicity m if and only if the polynomial 



(41) 



;=0 
passes with multiplicity m through the point (a, 7) = (pg{(x., /5), 
where q)g is the birational mapping defined in (|7]i. 

Proof. In order to prove the proposition, it is necessary and 
sufficient to show that 

F«,/3(S'(x,y)) = ii,,^{n{x,Y)) 

Note that in view of (ITTT i and (|38] l. we have ^{oc) 7^ when- 
ever g{a) ^ 0. Hence }ina{\p{X)') = 0. Let us define 



W(X,Y) 



!/.(X) 



Then^«,^(Q'(X,y)) = ii^^fi{U{X,Y)) in view of ©. There- 
fore, it is enough to show that 

}i,,l,{U{X,Y)) = ii,,^{n{X,Y)) 

But this follows from Theorem[T] since ZY(X, Y) does not have 
a pole at (a,^) and-H(X,Y) = 0g{U{X,Y)) in view of gD 
and (l42l) . where <P^ is the map introduced in Definition^ | 

Proposition|6]will be useful for those points (a, )S) £ V' for 
which g{a) 7^ 0. However, V' may also contain points whose 
X-coordinate coincides with the X-coordinate x, of one of the 
points in the re-encoding set TZ. For such points (a,/S) € V, we 
have g{oc) = by ( |37] ). In this case, we will use Proposition]?] 

Proposition 7. Let the point (a, /5) e "P' be such that a = x, for 
some point [xi,yi) in the re-encoding point set TZ. Then the 
polynomial Q' (X, Y), as defined in ( |40] |. passes through {a, /5) 
with multiplicity m if and only if the polynomial 

U'{X,Y) t^ Y. ^;(X) (X-a)^'-^' f;(X) Yi (43) 

passes with multiplicity m through the point (a, 7), with 7 given 
by7 = [i/g'{a), whereg'{X) is the derivative of g{X) in ( l37b . 

Proof. Again, to prove the proposition, it is necessary and 
sufficient to show that 

F«,/3(S'(X,Y)) = ^,,^{'H'{X,Y)) 

In order to do so, let us first introduce the auxiliary polynomi- 
als h{X) and (p{X) defined as follows: 



,(^^def |P0 ^ J^(x-X,) 



^(X) 



X-a 



def ^p{X) 






(X - a.yi 



n(x-x,y 



(44) 



(45) 



1=1 

Xj ^a 



Observe that, since a is a simple root of ^(X) by (|37| i, we have 
h{a.) = g'{oc), where ^■'(X) is the (first-order Hasse) deriva- 



tive of g{X). Hence 7 can be written as fi/h{a) and we have 
(a, 7) = cp],{ci,fi), where ip/, is the birational mapping in O 
defined with respect to h{X) rather than g{X). Also observe 
that (p{ci) 7^ in view of ( |45] |. Now, let us rewrite the expres- 
sion gii for Q'(X,Y) as follows: 

Q'{X,Y) = 4,{X)t'^j{X){X-ay'-hjiX) ^^ 

The rest of the proof is exactly the same as that of Proposi- 
tion|6l We define the rational function 



V(X,Y) 



dj Q'{X,Y) 

4,{X) 



and observe that f/fl,fi( Q'(X,Y)) = fi^fl(V(X,Y)) in view of 
dSll and the fact that }ia,ii (</'(X)) = 0. Further, V(X, Y) doesn't 
haveapoleat (a,^)and'H'(X,Y) = $;,(V(X,Y)), where ^;, 
is the map in (|9]l defined with respect to h{X). The proposition 
now follows immediately from Theorem[T] | 

Observe that the polynomials n{X,Y) and 'H'{X,Y) defi- 
ned in fiH i and ( l43T l are related via the transformation 



n'{x,Y) = {x-ocy'u{x, 



x-a 



Therefore 'H'{X,Y) passes through a given point (a, 7) with 
multiplicity m if and only if 

(46) 



H„,^(^(X-«)^'^(x,3^)j ^ m 

It follows from Corollary|5] Proposition|6] Proposition]?] and 
( |46] l that all the interpolation constraints of the shifted interpo- 
lation problem IPii^_i{V', M) can be expressed in terms of 
constraints on the polynomial 'H{X,Y) in dTTT l. 

It will be convenient to recast these constraints on 'H{X,Y) 
in terms of a point set V* and a corresponding set of multipli- 
cities M* . Given the original point set V in ( fTSi ) and the re- 
encoding point set TZ in ( i23T l, let us define 

def 
S = the set of points in T^\TZ whose X-coordinates 

differ from those of the re-encoding points 

def 
T = the set of points in T^\TZ whose X-coordinates 

coincide with those of the re-encoding points 

Thus the sets S and T form a partition of T^\TZ. The coordi- 
nate transformation consists of converting 'P\7^ into the redu- 
ced point set V* = {{xi,Zi) : {xi,yi) E P\7^}, where 

Vi - e{xi) 
def 



g'i^i) 



if {xi,yi)eS 



(47) 



if (xi,yi)eT 



We let S* and T* denote the sets of points in P* transformed 
from the points in S and T, respectively. Further, let M* C M 
denote the multiplicities of the points in P*, which are equal to 
the multiplicities of the original points in T^\TZ. Thus 

Propositions]6]and]7] and the coordinate transformation ( l47] i, 
are key to converting the original interpolation problem into the 
reduced interpolation problem, defined in the next subsection. 
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C. Reduced Interpolation Problem 

Before introducing the reduced interpolation problem in Defini- 
tion|6] we need to establish one more result. 

Lemma 8. Let Q'{X, Y) be a solution to the interpolation prob- 
lem lPii,_i{V', M). Letip{X) and'H{X,Y) be the polynomials 
defined in ( [38] l and (l4Tl i. respectively. Then 

degi,_i Q'(X,Y) = deg!/;(X) + de^^_{H{X,Y) 

Proof. Reformulating ( |40] l once again, let us write Q' (X, Y) 
as Q'(X,Y) = E-=oP;(^)^^'' where 



def ^{X)qj{X)tj{X) 



for 7 = 0, 1, . . . , r 



Note that Pj{X) is a polynomial (rather than a rational func- 
tion) for all /, and its degree is given by 

degpjiX) = degipiX) + degqjiX)tj{X) - jk 

since deg^(X) = fc by (l37T i. Since the r terms in E;=o Pii-^)^^ 
have distinct Y-degrees, it follows that 

degi^,_i Q'(X,Y) = 

= max.{degP;(X)+;(/c-l)} 

= degipiX) + max{deg^^(X)f^(X) -;} 

= degt/.(X)+degi_i?^(X,Y) | 

Lemma[8]is the reason we use the (1, — l)-weighted degree, 
rather than (1, fc— l)-weighted degree, in the definition below. 

Definition 6 (Reduced interpolation problem). Given the set 

ofpointsV* = {(xi,Zi), (x2/Z2),- ■ ■, {Xs^^rZ^^i^)}, partition- 
ed into subsets S* and T*, and the corresponding set of multipli- 
cities M* = {nix^,zi, nix2,z2i ■ ■ ■ I ^Xg_i^,Zs_i^}, ths reduced inter- 
polation problem consists of computing a nonzero polynomial 
T-L (X, Y) which can be expressed in the form 



n{X,Y) = Y.qj{X)tj{X)Yi 
;=0 



(48) 



and satisfies 

pixi,zi{n{X,Y))^m,^,zi V(x,-,z,)e5* (49) 

li,J{X-Xiy-'H(^X,^l^)^mx„zi V(x„z,)er* (50) 

such that degj _tT-L{X, Y) is minimal among all bivariate poly- 
nomials that satisfy the constraints (081), ( |49l ), and dSOl l. 

We shall refer to the reduced interpolation problem above as 
RIPj _l(P*, M*). The following theorem summarizes our re- 
sults and establishes the connection between RIPj _i(P*, M*) 
and the original interpolation problem IPj j;_i('P, M). 

Theorem 9. Let 7i{X,Y) be a solution to the reduced interpo- 
lation problem RIFj _i(P*,M*). Then a solution to the origi- 
nal interpolation problem IPi j--i i'P, M) is given by 



Q{X,Y) = ipiX)'H{X, 



giX) J 



(51) 



, def 



Proof. Let Q'(X, Y) = ip{X)'H (X, Y/g{X)). It is enough 
to show that Q'{X,Y) is a solution to the shifted interpolation 



problem IPj j;_i (P', M), since then the claim follows by Theo- 
rem[3] Proposition|6]and Proposition!?] along with ( |49] l and (|50] l 
imply that Q'{X,Y) passes through the points (xj.+i,i/i:+i), 
{^k+2>yk+2)> • ■ •/ (^s/i/s) in 'P' with the prescribed multiplic- 
ities. Using the definitions of ^(X), i/'(X), and t/(X) in ( |37] |. 
\, and ( [39] l, we can express Q'{X,Y) as follows; 



Q'{X,Y) = ^{X)'H(x, 



g{X) 



l^^cjjixMx)tj{x)^-^^ 



Y V 



(52) 



(53) 



= E^iiX)Ui^-^'t-~'^"^' (54) 

where the last equality follows from the fact that, for all inte- 
gers a and b, we have a — b= [a — b]^ — [b — a]^. It follows 
from ( |54] i and Corollary|5] that Q'{X,Y) also passes through 
the points {xi, 0), (x2, 0), . . . , (x^-, 0) in V' with multiplicities 
vi, V2, . . . , v/,. Thus Q'(X, Y) is in the ideal I{V', M) of poly- 
nomials that satisfy all the constraints of IPj /f_i(P', M). 

It remains to prove that S'(X, Y) is of minimal (1, fc— 1)- 
weighted degree in this ideal. Assume to the contrary that there 
exists a polynomial A{X,Y) E X{V',M) such that 

degi ,_i A{X, Y) < degi^,_i Q'(X, Y) (55) 

Combining Corollary|5]with the definitions of ^(X), i/'(X), and 
tj{X), we conclude that A{X,Y) can be expressed as 



A{X,Y) = ,/,(X)^a^(X)t^(X)(^-^^ 



(56) 



Given ([561, we can construct a solution to RIPj _i {V*, M*) as 
follows. Define B{X,Y) = A{X,Yg{X))/ip(x), namely 



B{X,Y) = ^fl^(X)f/X)Y/ 
;=0 



(57) 



Note that ( |57] ) is of the form ( |48l ). Furthermore, Propositions |6] 
and|7]imply that 6(X, Y) satisfies (gill and (IIOI). By Lemma[8] 

degi,_i^(X,Y) = deg!/^(X) +degi^_ifi(X,Y) 

which, in conjunction with ( fSSl ), implies that deg^ _2S(X, Y) 
is strictly smaller than degj _i'H{X,Y), a contradiction. | 

The reduced interpolation problem RIPi _i {V*, M*) can be 
solved using a variant of Koetter's interpolation algorithm, de- 
scribed in the previous section. We need to make the following 
three modifications. First, in view of (l48T l. we need to initialize 
the Grobner-basis polynomials as follows: 



gp{X,Y) := tj{X)Y' 



1=1 



,\j~^i]' 



(58) 



for j = 0,1, . . .,r. This initialization guarantees that the Grob- 
ner-basis polynomials satisfy (l48l l through all the iterations of 
the algorithm. Second, in view of ( |50l ), in those iterations that 
deal with points (a, 7) G 7"*, the discrepancies A, need to be 
defined differently. Specifically, suppose that a = x, for some 
point (x,, y,) in the re-encoding point set TZ, and let v, denote 
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the multiplicity of this point (as before). Then the discrepan- 
cies should be computed as follows: 



A; 



coef X^'g,- X + a, 



y + 7 

X 



;X''yf 



(59) 



for ;' = 0,1, ...,r. Third, instead of the monomial order -</( 
defined in ( |20] i. we need to order the bivariate monomials as 
follows: we say that X''Y^-<_i X^Y^ iff 

(fl — b < /— ;') or {a — h = i — j and b < j) (60) 

Note that -< _i is not a monomial order, since it does not have 
the well-ordering property (there is no smallest monomial un- 
der -<-i). This, however, has no effect on Koetter's algorithm, 
and we may safely replace -<k with ^_i throughout. 

Example 2c. Consider again the re-encoding transformation in 
Example|2al Recall that k = 2, and the re-encoding set TZ con- 
sists of (xi,i/i) = («,«"*) and (x2,i/2) = (a^,a^), with multi- 
plicities 1/1=2 and 1/2 = 1 respectively. Thus 

2 



g{X):= UiX-^') = (X-<^)(X-a2) 



(61) 



i=l 



ip{X) 



x,y 



(X-a)2(X-a2) (62) 



according to ( |37] l and ( |38] l, respectively. Since r = 3, there are 
four tail polynomials fo(X), ti(X),f2(X),f3(X). However, the 
first two are trivial: tQ{X) = fi(X) = 1. Indeed, let us define 
i^min = iriini<j,<jj- Vj. Then the definition of tail polynomials in 
(|39] l implies that f;(X) = 1 for / ^ v^j^. In fact, if Vmm ^ ''. 
then all the tail polynomials are trivial, and the constraint (|48T l 
on the form of a solution to RIPi _i(P*, M*) becomes vacu- 
ous. In this example, however, we have Vmin = 1^ and 

X-a^ 



(X-a)(X-a2)2 = oc^ + a^x + aX^ + X^ 



h{X) 

f3(X) 

Thus, according to (ISST l. the Grobner-basis polynomials should 
be initialized as follows 

'(0)/ 



^o(X'^) 



q"'(X,Y) 
G2\X,Y) 

g^yx,Y) 



fo(X) = 1 

h{x)Y = y 
f2(x)y2 = (X- 



a2)y2 

a^X + aX^ + X^)Y^ 



f2(X)y3 = (, 

Next, let us compute the reduced point set V* using the coor- 
dinate transformation defined in ( |47] i. The partition of the five 
points in V\TZ into the sets S and 7" is given by 



S 
T 



{(«3,l),(a3,a),(l,a),(l,l)} 
{(<.l<.3)} 



since ^(a^) = 0, while g{a.^) = oc^ and g{l) = cp- are nonzero. 
Also g'(a^) = OL^- Dividing the shifted interpolation points in 
Examplel2bl bv a^, a^, or a^, as appropriate, we compute 






{(a3,a2),(a3,a3), (1,0), (!,«)} 
{(«M)} 



The multiplicity of all the points in V* = S* \J T* is 1, as in 
Examplel2bl This determines all the interpolation constraints of 



RIPl^_l(P*,M*), and Koetter's interpolation algorithm, sub- 
ject to the modifications (|59] l and (|60] l, now proceeds through 
the 5 iterations detailed in Table III. Again, at each iteration, 
the Grobner basis polynomials are arranged in Table III in as- 
cending order with respect to ^_i. Thus the output of the al- 
gorithm is the polynomial 'H{X,Y) = t/^ (X, y), given by 



n{X,Y) = {ct' + X)Y + (a5 + a^X + x2)y2 



(63) 



It is easy to verify that ip{X)'H (X,Y / g{X)) , where g{X) and 
ip{X) are given by ( l6Tl i and ( |62] |. is precisely the polynomial 
Q'(X, y) in ( |33] |. It can be also verified directly that 

is precisely the solution to the original interpolation problem 
found in Example|2al (cf. ( [3T| i or the last row of Table|I]i. □ 

Remark. Strictly speaking, the set J{V*, M*) of all bivari- 
ate polynomials over Wq that satisfy the constraints ( |48] ), ( |49l ). 
(ISOl l is no longer an ideal of F^ [X, y] . The problem is with the 
constraint (gS)!. For example, if A{X,Y) e J{V*,M*) then 
y.4(X, y) is not necessarily in J{V*, M*) since it may not be 
possible to express y.4(X, y) in the form of ( |48] l. Neverthe- 
less, if we take the intersection of J {V* , M*) with the set of 
polynomials in F^ [X, y] whose y-degree is at most r, then this 
intersection is a module over F,j [X] . What the (modified) Koet- 
ter algorithm computes is the Grobner basis for this module. 

In summary, we emphasize that although RIPj _i(P*, M*) 
appears to be much more convoluted than the original prob- 
lem Wif^_i{V,M), its complexity is often orders of magni- 
tude lower. This is due to the fact that we do not even need to 
consider the first k points of V in computing Ti. (X, y). In other 
words, these k interpolation points (which can be selected to 
have the largest multiplicities) effectively disappear. The com- 
plexity savings are not apparent in the toy example studied in 
this section. However, they become prominent for long high- 
rate Reed-Solomon codes of practical interest (cf. Example I). 



V. The Factorization Procedure 

The reduction in complexity achieved in the previous section 
would be much less significant if one had to actually use ( BTI ) 
to compute a solution to the original interpolation problem (or 
the shifted interpolation problem), and then factor this solu- 
tion in order to recover the transmitted codeword. The problem 
is that ( BTI ) involves multiplication by ip{X), which is a large 
polynomial, usually much larger than the solution T-i{X, Y) to 
RIPl _l(P*,M*) computed by the Koetter algorithm. For in- 
stance, in the situation of Example 1, typical for Reed-Solomon 
codes of practical interest, we have degi/'(X) = 1663 whereas 
degi o'H(X, y) = 44 and degg ^^^(X, y) = 6. Thus in factor- 
ing "H (X, y), one deals with polynomials of degree at most 44, 
whereas factoring Q'(X,y) = ip{X)'H{X,Y/g{X)) involves 
processing polynomials of degree 1598. 

Fortunately, the following theorem shows that instead of fac- 
toring Q(X, y) or Q'(X, y), we can directly factor the much 
smaller polynomial ?^ (X, y ) to find the transmitted codeword. 
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Theorem 10. Let 'H{X,Y) be a solution to the reduced inter- 
polation problem RIPi _i(P*,M*) and let 

y-g(x)\ 



Q{X,Y) = ip{X)n{X, 



(64) 



Ycr{X)-co{X) 



(7(X) Y 



giX) J 

be the corresponding solution to the original interpolation prob- 
lem IPi,jc_i {V, M). IfQ{X, Y) has a factor Y - f{X), where 
deg/(X) < k, then 'H{X,Y) has a factor of the form 

where a{X) andco{X) are the error-locator and erwr-evaluator 
polynomials for the re-encoding positions with respect to f{X). 

Proof. Recall that the re-encoding point set TZ consists of 
the k points {xi,y\), {x2,y2),- ■ ■, {^kiyk) with Xi,X2/- ■ ■ ,^k 
being distinct elements of D. By the error-locator polynomial 
for the re-encoding positions with respect to /(X), we mean 
a monic polynomial (j{X) with cr(x,) = whenever /(x,) 7^ 1// 
and a{x) 7^ for all other x G Fg. In other words, let 



£ ^^ {ie{l,2,...,k} : y,7^/(x,)} 



(66) 



be the set of error locations among the re-encoding positions. 
Then the error-locator polynomial (t{X) is given by 

a[X)^MYl{X-x,) (67) 

Further, let S'^ be the complement of £ in the set {1, 2, ... , k}, 
and define the polynomial 

X{X)''M Y[{X-Xi) (68) 

By an error-evaluator polynomial for the re-encoding positi- 
ons with respect to /(X), we mean a polynomial Cjo{X) such 
that whenever /(x,) 7^ i/,-, their difference is given by 



f{^i)-yi = a;(x,-)A(x/) = w{xi) 



(j'{Xi) 



(69) 



where g'(X) and cr'{X) are the (first-order Hasse) derivatives 
of the polynomials ^(X) in dSTJ l and (t{X) in ( |67] |. 

Now suppose that Q(X, Y) has Y — /(X) as a factor. That is, 
there is a polynomial A{X,Y) such that 



S(X,Y) = Y-/(X)L4(X,Y) 



(70) 



Let e(X) be the re-encoding polynomial defined in (|24| ). and let 
Q'{X,Y) = Q(X,Y + e(X)) be the con-esponding solution to 
the shifted interpolation problem. It follows from (ITOl i that 

Q'{X,Y) = (Y-l^{X)')A{X,Y + e{X)) (71) 

where t]{X) =^ f{X) - e{X). If n{X,Y) is the con-esponding 
solution to the reduced interpolation problem, then 



^p{x)n(x, 



Y 



g{X) 



Y->]{X))A{X,Y + e{X)) (72) 

Finally, applying the birational transformation <Pg in (|9| to both 
sides of ( |72| |. we obtain 

xpiX)'H{X,Y) = (Yg(X)-;/(X))fi(X,Y) (73) 

where B{X,Y) = A(X,Yg{X) + e{X)). It follows that the 
polynomial ip{X)'H{X,Y) is divisible by Yg{X) - t]{X). 



Now observe that T]{xi) = /(x,) — e(x,) = /(x/) — y; = 
for all / G £'^. Therefore ^(X) is divisible by A(X), and we can 
express i]{X) as co{X)A{X) for some polynomial ^^(X). But 
then 7^ (X;) = f{xj) — y,- = a;(x/)A(Xj) for all z, and therefore 
a;(X) is an error-evaluator polynomial with respect to /(X), 
as defined in i6% . Furthermore, since ^(X) = cr{X)A{X) and 
t]{X) = co{X)\{X), we can rewrite (ITsT l as follows: 

ip{X)n{X,Y) = A(X)(Y(7(X)-a;(X))s(X,Y) (74) 



Notice that //(x,) 7^ for all i e £, and therefore a;(x,) 7^ 
for all i E £ as well. Referring to ( |67] ). this implies that (j(X) 
and co{X) are relatively prime. Consequently, the polynomials 
!/?(X) and Ycr(X) — co{X) are also relatively prime. With this, 
it now follows from ( TMb that Ycr(X) — cv{X) must be a factor 
of'HiX,Y), as claimed. | 

We shall see shortly how the polynomials cr(X) and co{X) 
can be exti-acted from a solution V. (X, Y) to RIPi^_i (P*, M*). 
For the time being, let us assume that these polynomial are al- 
ready known, and explain how decoding should be completed. 
First, we find the roots of (j(X), using either Chien search [15, 
p. 276] or direct factorization [21]. This reveals the set of error 
locations £ in ( |66] |. Next, we use a;(X) along with the deriva- 
tives of cr{X) and g{X) to compute the error values 

g'i^i) 



def 



co{xi) 



a'{x,) 



for all i e £ 



(75) 



Since the set 72. = {{xi,yi),{x2,y2), ■ ■ ■ , (^Jcl/J:)} is known to 
the decoder, this makes it possible to determine the values of 
the polynomial /(X) at k distinct points, namely 



/(^O 




(76) 



Since deg/(X) < k, the entire polynomial /(X) can be recov- 
ered from these k values using the same procedure (e.g., Lag- 
range interpolation) that was used to compute e(X) in ( |24] |. 

It remains to show how to compute the polynomials cr{X) 
and iv{X) from 'H{X,Y). Note that the ratio a;(X)/tr(X) is 
a rational function, which can be expanded as the power series 

a;(X) def 



To + liX + 12X^ 



Et.x^ 



(77) 



provided cr(0) 7^ 0, as we assumqj Roth and Ruckenstein [18] 
view a bivariate polynomial A{X, Y) as a polynomial in Y with 
coefficients in the ring Wq [X], and define a Y-root of A{X, Y) 
as any element /(X) of Fg[X] such that ^(X,/(X)) is iden- 
tically zero. They develop an efficient iterative algorithm for 
computing the Y-roots of a given bivariate polynomial. The 
Roth-Ruckenstein algorithm of [18] takes d iterations to pro- 
duce all the Y-roots of degree at most d. Herein, we observe 
that the Roth-Ruckenstein machinery [18] extends to rational 
functions. Namely, let us extend the definition of a Y-root of 

' We observe that i7(0) = if and only if Xj = for some re-encoding point 
{xj, yi) with i £ £. This situation can be avoided by not including points of type 
(0,1/) in the re-encoding point set 7?. (in fact, this will happen automatically if 
the code support set P is a subset of F,,"). If one insists on including a (0, y) 
point in 7?., we can expand w{X)/a(X) in a power series about another point 
a with 17(a) 7^ 0. In this case, a method similar to the one we describe applies. 
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A{X, y) to mean any element 7(X) of the field Fq(X) of rat- 
ional functions over F^, such that A[X,j{X)^ is identically 
zero. Then d iterations of the Roth-Ruckenstein algorithm pro- 
duce the first d coefficients of the power-series expansion of 
any such Y-root. Since a;(X)/cr(X) is a Y-root of 'H{X,Y) 
by TheoremfTol we can use the Roth-Ruckenstein algorithm to 
iteratively generate the coefficients 70, 7i, 72, ... in dTTb . 

How many coefficients do we need? This depends on the 
degree of c{X), which is equal to the number of errors in the 
re-encoding positions. Note that cr{X) and a;(X) can be found 
from the power-series expansion of co{X) /cr{X) by a Pade ap- 
proximation procedure, such as the Berlekamp-Massey algo- 
rithm. It is well known [15, p. 366] that at most 2f coefficients 
of the power series are required to reconstruct <jo{X) and cr{X), 
where t = deg(j(X). In theory, degcr{X) ^ fc in view of (l66ll 
and i&U . so 2k iterations of the Roth-Ruckenstein algorithm 
always suffice. In practice, we would expect much less than k 
errors in the k re-encoding positions. Indeed, in soft-decision 
decoding, we usually choose the re-encoding point set to con- 
sist of the k most reliable positions. In this case, the number of 
errors in these positions will be small with high probability, 
since most errors occur in the n — k least reliable positions. 
Consequently, in practice, one would run the Roth-Ruckenstein 
algorithm for 2t iterations, where the integer parameter t ^ A: 
is a pre-determined bound on the number of errors we expect 
(with high probability) in the k re-encoding positions. 

Our results in this section can be summarized in the form 
of a reduced factorization algorithm, presented below. 

Reduced Factorization Algorithm 

Input: A solution 'H{X, Y) to the reduced interpolation prob- 
lem RIPi _i {V*, M*), an integer parameter t, and the re- 
encoding point set Tl = {{xi, 1/1), {x2, yi),---, {xk, i/jt)} ■ 

Factorization: Apply the first 2t iterations of the Roth-Ruck- 
enstein algorithm [18] to ?:^(X, Y), thereby producing the 
first 2t coefficients 70, 71, . . . , 72t-1 in the power-series 
expansion of the Y-roots of 'H{X,Y). There are at most 
r such Y-roots, where r = degg^?^(X, Y). 

Fade Approximation: Treat each sequence 70, 71, . . . , 72t-i 
produced in the previous step as a "syndrome" sequence; 
then compute the coefficients ctq, a^, . . .,at of the shortest 
linear feedback shift-register that generates this syndrome 
sequence (Berlekamp-Massey algorithm). Next, find the 
coefficients coq, coi, . . . ,cot of the error-evaluator polyno- 
mial a;(X), e.g., using the convolution: 



CO; 



^o7; + o"i7!-i 



^i-l7l + ^/7o (78) 



Root Finding: Find all the roots of the error-locator polyno- 
mial cr(X) = (7o + (JiX + ■ ■ ■ + a'tX^ and for each such 
root compute the corresponding error value using dTSl ). 

Corrected Re-Encoding: Set the k corrected re-encoding val- 
ues /(xi),/(x2)/ ■ • ■ rfi^k) as in ( l76l ). and find the unique 
polynomial /(X) of degree < k that agrees with these val- 
ues (using, for example, Lagrange interpolation). 

Output: Return the list of all such polynomials /(X). 



Not all the Y-roots produced by the 2t iterations of the Roth- 
Ruckenstein algorithm necessarily correspond to a valid cr(X) 
and co{X) pair Just as in conventional Berlekamp-Massey de- 
coding, we expect the error-locator polynomial (j(X) and the 
error-evaluator polynomial co{X) to satisfy certain conditions. 
If these conditions are violated, the corresponding Y-root can 
be rejected. For example, a Y-root sequence 70, 71, . . . , 72t-1 
can be safely rejected if any of the following occurs: 

a) the degree t of cr{X) is strictly greater than t; 

b) the convolution (fTSl l produces nonzero values for / > t; 

c) the polynomial cr{X) has less than t distinct roots in F^; 

d) an error value e,- = is computed according to dTSl l. 
We leave the proof of the underlying properties of (t{X) and 
a;(X) as an exercise for the reader. In practice, all these con- 
ditions are used to detect false Y-roots in the reduced factor- 
ization algorithm. As a result, with very high probability, only 
one polynomial /(X) is returned by the algorithm. 

Example 2d. Consider again the situation in Example 2. Recall 
that the solution Q{X,Y) to the original interpolation prob- 
lem, found by the Koetter algorithm, factors as follows: 

S(X,Y) = {a.^+X){Y-ia^+a^X)){Y-{a^ + ci^X)) 

The re-encoding point set in Examplel2al consists of the points 
(xi,i/l) = (a, a'^) and (^2,1/2) = (a^, a^), and the resulting re- 
encoding polynomial is e(X) = a^ + a^X. Solving the reduc- 
ed interpolation problem produces the polynomial 

'H{X,Y) = {a.^ + X)Y + {oc^ + a^X + X^)Y^ 



as shown in Example|2c] (cf . (I63l l and the last row of Table III). 
Observe that this polynomial factors as follows: 

n (X, Y) = (a^ + a^X) Y ( Y(l + a^X) - a^) (79) 

Applying the Roth-Ruckenstein algorithm [18] to ?^(X,Y) pro- 
duces two syndrome sequences. The first sequence, which cor- 
responds to the factor Y in (l79T l. is identically zero. The short- 
est feedback shift-register that generates the all-zero sequence 
is (j(X) = 1, and the corresponding error-evaluator polynomi- 
al is co{X) = 0. This indicates the absence of errors in the re- 
encoding positions, and the resulting message polynomial is 



/i(X) 



eX) 



■ a°X. Note that Y - /i(X) is, indeed. 



a factor of Q{X,Y). The second syndrome sequence, corres- 
ponding to the factor Y(l + a X) — a. in ( l79l l. is given by 

70/ 71/ 72/ • ■ • = "i^, a''/ '^z cic^, Dc^, 0?, 1/ oc" ,... (80) 

Applying the Berlekamp-Massey algorithm to the syndrome se- 
quence dSOl l. we find that it is generated by cr(X) = 1 + oc'X. 
The convolution dTST l then produces a;o = oc" and <jO\ = 0, so 
that a;(X) = ct'. The error-locator polynomial has a single root 
at X2 = i^- Thus the set of error locations is £^ = {2}, and the 
single error value is given by 



£2 



a;(x2) 



r'(^2) 



In accordance with ( l76l ). we now set JTJyXj) = yi^ ei = 0? 
and f^ixx) =y\= oc . The unique polynomial of degree < 2 
that has these values is fiiX) = a^ + iX^X. This polynomial, 
indeed, corresponds to the other factor of Q{X, Y), and to the 
transmitted codeword (cf. Example|2a]i. □ 



12 



Submitted to the IEEE TRANSACTIONS ON INFORMATION Theory, April 2010 



VI. Conclusions 

A complete proof of the applicability of the re-encoding and 
coordinate transformation method to the bivariate polynomial 
interpolation process in algebraic list-decoding of Reed-Solo- 
mon codes is presented. A detailed example is given to illust- 
rate the entire re-encoding and coordinate transformation pro- 
cess. In addition, it is shown how the factorization procedure 
can be modified to accommodate the reduced interpolation pro- 
blem. Thereby, factorization complexity is also decreased sig- 
nificantly since the required number of iterations of the Roth- 
Ruckenstein factorization algorithm is reduced from k to 2t, 
where t is a small number (we found that, in practice, t ^ 6 
often suffices, even for long high-rate Reed-Solomon codes). 
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(^,y) 


ntx.y 


Grobner Basis Polynomials 


(a, a") 


2 


go{x,Y) 

Qi(X,Y) 
Q2(X,Y) 
GiiX.Y) 


= a + X 
= a* + Y 

= ftS + Y3 






Qi(X,Y) 
Qo(X,Y) 
Q2(X,Y) 
Q3(X,Y) 


= a-* + Y 
= a2 + X^ 

= Ci^ + Y^ 






Qo(X,Y) 
Qi{X,Y) 
Q2(X,Y) 
Q3(X,Y) 


= a2 + X2 

= (a5 + a''X) + (a + X)Y 

= a + Y2 

= aY + Y'^ 






{a^,cc^} 


1 


Qi(X,Y) 
Q2{X,Y) 
Qo[X,Y) 
QiiX.Y) 


= (a*- + ot^X + ofiX^) + (ft + X) Y 
= (a3+a5x2)+y2 

= a* + K^X + a^x^ + X3 
= (a'^ + a^X^) +uY + Y^ 






(a2,a3) 


1 


Q2(X,Y) 
Qo(X,Y) 
Qi(X,Y) 
Qi(X,Y) 


= (a* -(- a'lX + aX2) + (a + X) Y + Y^ 
= a* + a2X + «2x2 + x3 
= (a + a2x2 + a«^x3) + (fl;3 + a^X + X2)Y 
= {o? + a^X) + (a^ + a5x)y + y3 






(a^l) 


1 


Qo{X,Y) 
Qi(X,Y) 
Q2{X,Y) 
Qi[X,Y) 


= (1 + a^X + X3) + («2 + aX)Y + aY^ 
= (a5 + fl;<^X + (x^X^ + a^X^) + (aX + X^)Y + a^Y^ 
= (1 + a^x + aX3) + (a^ + X + X2)Y + (a^ + x)Y^ 
= (a^ + a^X) + (ftS + a5x)Y + Y^ 






i'^',^) 


1 


Qi{X,Y) 
Q2{X,Y) 
Qi(X,Y) 
Qo{X,Y) 


= (a* + a^X + a^X^ + a^X^) + (a^ + X^) Y + a^Y^ 

= (1 + a^X + aX3) + (a* + X + X2)Y + (a^ + xjY^ 

= (a* -h a'-X^) + (a* + a^X) Y + Y^ + Y^ 

= (ft;3 -h a^X + a^x^ + a^X^ + X") + (a^ + aX + aX^)Y + (r 


* + aX)Y2 




(l,a) 


1 


Q2(X,Y] 
Q3{X,Y) 
Qo[X,Y) 
Qi(X,Y) 


= (1 + a^X -h aX3) + (a^ + X + X2) Y + (a^ + X)Y2 

= (a^ + X + aX2 -t- aX3) + (a + a^X + o?X^)Y + Y^ 

= (a2 + K^x + a'lX^ + X*) + {a2 + aX) Y + (1 + aX) Y^ 

= (a* + X2 + a3x3 + a^X*) + (a^ + ^Zx + X^ + X3)Y + (a* 


+ a4X)Y2 




(1,1) 


1 


Q2(X,Y) 
Q3{X,Y) 
Qi{X,Y) 
Go{X,Y) 


= (1 + a^X + aX3) + (a" + X + X2)Y + (a^ + X)Y^ 

= (a5 + X + aX2 + aX3) + (a + ft:*X + a^X^) Y + Y^ 

= (a* + X^ + a3x3 + a^x") + (a2 + aZ^ + x^ + X^jY + (a* + a^X) Y^ 

= (a2 + a^X + a<5x2 + a^X^ + X* + X^) + (a^ + a^X + OiX^)Y + (1 + a^X + 


aX2)Y2 



TABLE I 
Iterations of the Koetter Algorithm for the Original Interpolation Problem in ExAMPLE[2al 
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(^,y) 


mx,y 


Grobner Basis Polynomials 


(«,0) 


2 


£?o(x,y) = « + x 
£?i(x,y) = y 
&(xy) = y2 
fo(x,y) = y3 


go(x,y) = a + x 
gi(x,y) = (a + x)y 
&(x,y) = y2 
03(x,y) = y3 


6;o(x,y) = a2 + x2 
gi(xy) = (« + x)y 
&(x,y) = y2 
fo(x,y) = y3 


(0^,0) 


1 


eo(X, y) = a* + a^X + OL^X^ + X^ 

gi(x,y) = (a + x)y 
&(xy) = y2 
e?3(x,y) = y3 


{a^,a) 


1 


g2(x,y) = (a2 + ax)y + y2 

6?o(X, y) = (a^ + a^X + a2x2 + X3) + (a5 + a*X)y 

ei(x,y) = (a« + x + x2)y 

fo(X,y) = (a3+a2x)y + y3 


(^3,1) 


1 


eo(X,y) = {«^+a^X + a^X^ + X^) + [oi + X)Y + Ci^Y'^ 

gi(x,y) = (a'' + x + x2)y 

&(X,y) = (a5+aX + fl;X2)y+(a3 + X)y2 

e3(x,y) = {a2 + ax)y + a3y2 + y3 


(1,0) 


1 


gi(x,y) = (a'' + x + x2)y 

&(X,y) = {a5+aX + aX2)y+(a3 + X)y2 

e3(x,y) = (a2 + ft;X)y + a3y2 + y3 

go(X,y) = (a^+aX + a<'X3 + X*) + (a + a3x + x2)y+(a* + a''X)y2 


(1,«3) 


1 


&(x, y) = (a* + X + x2)y + {a3 + x)y2 

g3(x,y) = (a + a3x + x2)y + a3y2 + y3 

eo(X,y) = {a*+ttX + a<'X3 + X*) + (a + a3x + x2)y+(a* + a''X)y2 

ei(x,y) = (a* + a5x + x3)y 


(«2,<.4) 


1 


&(x, y) = (a* + X + x2)y + (oc' + x)y2 
Gsix, y) = (a + a3x + x2)y + a3y2 + y3 

Gi (X, Y) = (a* + ttX + a''X3 + X*) + (a^ + cP-x + X^ + x3)y + (a* + a*X)Y'^ 
Go(X, Y) = (a'^ + a'^X + olX^ + ffiX3 + X* + X^) + {i>? + c&X + a^X^ + X3)y + (a* 


+ a3x + a4x2)y2 



TABLE II 

Iterations of the Koetter Algorithm for the Shifted Interpolation Problem in ExAMPLElIal 
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{x,z) 


mx,z 


Grobner Basis Polynomials 


{a^,cc') 


1 


e2(x,Y) = ffiy + (a2 + x)y2 

e?o(X,Y) = l + «''Y 

gi(X,Y) = (a3 + x)Y 

g3(X,Y) = ii?Y + (a5 + a''X + aX2 + X=')Y=' 


(«3,«2) 


1 


go(X, Y) = 1 + Y + (a6 + «*X)Y2 

e?i(X,Y) = (a3 + x)Y 

^(XY) = (a* + aX)Y+ (a? ^ Ci^X^X^)'Y'^ 

g3(X,Y) =aY + (a5 + a3x)Y2 + (a^ + a''X + aX^ + X^jY^ 


(1,0) 


1 


£?i(XY) = (a3 + x)Y 

e2(X,Y) = (a* + aX)Y + (a^ + a5x + x2)y2 

g3(X,Y) = aY + (a5 + fl:3x)Y2 + (a^ + a^X + aX^ + x3)y3 

go(X,Y) = (l + X) + (1 + X)Y+ (a^ + a3x + a*x2)Y2 


(1,«) 


1 


e2(X,Y) = (a3 + X)Y+ (a5 + a5x-fX2)Y2 

g3(X,Y) = (1 + X)Y + (a5 + a3x)Y2 -f (a^ + a*X + aX^ + x3)y3 
eo(X,Y) = (l + X) + (1 + X)Y+ (a'' + a3x + a*x2)Y2 
ei(X,Y) = (fl;3 + aX + X2)Y 


(«M) 


1 


e2(X,Y) = (a3 + x)Y+ (a5 + a5x + X2)Y2 

g3(X,Y) = (1 + X)Y + (a5 + a3x)Y2 + (a^ + a*X + aX^ + x3)y3 

01 (X, Y) = (1 + X) + (a + a^X + X2) Y + (a*' + a^X + a*X2)Y2 

eo(X,Y) = (a^ + a^x + X^) + {a2 + fl;6x + x2)Y + [ct + aX + c&X^ + c&X^)Y^ 



TABLE III 
Iterations of the Koetter Algorithm for the Reduced Interpolation Problem in Example[2c1 



